The Influence of Cylindrical Inclusions on the Stability of a Directionally Solidified 

Interface 
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It is well known that the presence of an inclusion 
in the melt near a solidifying front induces a local 
deformation in the latter provided that the melt's 
thermal conductance differs from that of the inclusion 
Q, 1^ 1^ Q|. This local interfacial deflection is caused 
by the modification of the thermal gradient in the melt 
near the particle. The long time evolution of this local 
disturbance of the interface has been investigated for 
the case of a sphericalparticle in a pure substance 0, 
and in a binary alloy Q- It is discovered that, provided 
that the particle-interface distance falls below a critical 
value, the induced perturbation grows and destabilizes 
the whole solid-liquid interface. This newly uncovered 
morphological instability is manifested only for some 
combination of the physical and processing parameters, 
and its onset is attributed to the reversal of the thermal 
gradient in the liquid gap between the particle and the 
interface |^. This instability, whose characteristic size is 
of the order of the particle's diameter, occurs at pulling 
speeds that are below the threshold for the onset of the 
Mullins-Sekerka instability. 

In reality, however, the inclusions have approximately 
cylindrical shapes. For example, in their study involving 
a particulate metal matrix composite of magnesium 
base allov that is reinforced with SiC particles, Essa 
et. al state that the inclusions resemble cylinders 
whose longitudinal sections can be approximated by 
rectangles of sides a and b {a ^ b). The pupose of this 
communication is to extend the stability analysis that 
was carried out for a spherical inclusion in Q to the 
more realistic case of a cylindrical inclusion. We let 
the inclusion have the more general shape of an elliptic 
cylinder, i.e. the longitudinal section is an ellipse of 
semi- axes a and b as illustrated in figure 1. The effect of 
the aspect ratio c = b/a, which emerges in the analysis 
as an important factor, will be examined. 

Consider a system in which an inclusion is immersed in 
a bath of pure melt that is undergoing directional solidi- 
fication. We ignore the gravitational effects and assume 
that the inclusion is both solid, i.e. nondeformable, and 
insoluble in the melt. Let the inclusion's position be such 
that the distance from its center to the planar solid-liquid 
interface is Zq, and its axis of symmetry is parallel to the 
y— axis. The problem, being invariant with respect to 
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FIG. 1: A schematic diagram of a vertical solidification setup 
of a pure melt in which is embedded a cylindrical inclusion 
of aspect ratio c = b/a. The distance from the inclusion's 
center to the planar interface is Zq. The cylinder's axis of 
symmetry is parallel to the V— axis, and the planar interface 
coincides with the X — axis. The front's growth rate is V, 
G'^' is the imposed thermal gradient, and Ni is the normal 
to the interface pointing into the melt. 



the y— direction, is thus two-dimensional. The physi- 
cal process is governed soleley by heat diffusion in the 
melt, solid, and inclusion with appropriate conditions at 
the interface, at the inclusion's surface, and at infinity. 
Let T'''^\X, Z) represent the temperature in phase (g), 
where the superscript {q) denotes symbolically the liquid 
phase, (g = L), the solid phase, {q = S), and inclusion, 
{q = I), and the couple {X, Z) represents the Cartesian 
coordinates of a point in the XZ-plane of a right-handed 
coordinate system relative to which the solid-liquid in- 
terface is at rest at Z = 0, i.e. the X-axis. The melt 
occupies the top region, Z > as shown schematically in 
Fig. 1. The equation of the cylindrical inclusion, whose 
center is at the point (0, Zq), is given by 
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where a and 6 are the semiaxes in the X and Y directions, 
respectively.. The nondimensionalization that is adopted 
here makes use of a, a/V, and the melting point r„i as 
scales for length, time, and temperature. The follow- 
ing dimensionless system of equations and corresponding 
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boundary conditions emerges, 



dt dz 



y(s) ^ ^ 1 _ CTK, at z = 0, 



(2) 
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Svn^k— , at z = 0, (4) 
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T(^)=T(^), JJ^^^lll, on 5/ (5) 



dnp dnp 



dz 



dz 
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The symbols that appear in Eqs.(2)-(6) are defined 
in the foUowing: A stands for the two-dimensional 
Laplacian, dxx + dzz, e = aV/D^^^ is the Peclet number, 
and A(9' = D'^'>^/D^^\ where ZJ^') ig the thermal 
diffusion coefficient in phase (q). Equation (3) describes 
the continuity of temperature at the interface, with an 
interface equilibrium temperature that accounts for the 
Gibbs-Thomson effect. The parameter a = asL/o-L 
is the surface energy parameter, with crg^ being the 
interface excess free energy, L is the latent heat of fusion 
per unit volume, and k is the front's curvature consid- 
ered to be positive when the sohd bulges into the melt. 
Equation (4) describes the heat balance at the interface, 
where S = LVa/{Tmk^^^) is the Stefan number, d„ is 
the normal interface velocity, k — k^^^ /k'^^\ and d/drii 
is the normal derivative at the solid-liquid interface. 
Equation (5) represents the continuity of temperature 
and of the heat flux on the inclusion's boundary dl, 
where a = k'^-^^ /fc^^' is the thermal conductance ratio, 
and d/dnp is the normal derivative at the inclusion's 
surface. Finally, Eq. (6) describes the far field con- 
ditions for the temperature, where Gs — aG'^"^-' /Tm 
and Gl = aG^^^T^, with G^^) and G^^) being the 
externally imposed (dimensional) thermal gradients in 
the liquid and solid phase, respectively. 

The average size of the reinforcements that are typi- 
cally used in the manufacture of PMMCs is about 12 m 
0, while the melt's thermal diffusion coefficient is of the 
order lO^^m^/s 0. Thus, if we assume that the growth 
rate is low enough, say V « 10~®m/s, then the corre- 
sponding Peclet number e w 10^^. Note that the small- 
ness of e is due to the fact that it is proportional to the 
product of the velocity and the particle's radius, both of 
which are considered very small. Thus, the steady state 
form of Eq. (2) satisfies 



dz 



(7) 



As in ref. we could take advantage of the small- 

nes of e and consider the expansions, with e <C 1, 
T(9) = Ts^''^ + e6i(9) -I- 0(e2) and, for the interface, 
z — + £?/ + O(e^). At order unity, the base state 
satisfies Laplace's equations with coupled boundary 
conditions at the inclusion's surface and at the planar 
interface. This problem is mathematically tractable, and 
the next order perturbations are easily calculated. We 
refer the interested reader to for details. However, 
for both simplicity in the calculations and clarity in the 
presentation, we set e = in the present work. This 
simplification will not alter the key qualitative features 
of the instability mechanism that is the subject of this 
communication. 

It is known that the thermal field in an infinite medium 
in which is embedded an inclusion of different material is 
given by [3, 



— — — (z + nHo)Cx(x,z) , 



where 



C 



X{x,z,n) 



du 



2 7a(...,„) (c2-fj.)3/2(l+^)l/2' 



(8) 
(9) 



c = b/ a, Hq = Zq/q, and X{x, z, n) is the positive root of 
the equation, 

+ iiJl^ ^ 1. (10) 



1 -|- X{x,z,n) (? -|- \(x,z.n) 

For an elliptic cylinder, Eq. (9) reduces to 



C 



1 -I- A(a;, z, n) 



while for a circular cylinder, we have 



1], (11) 
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A(a;,js,n) 



2[a;2 + (z - uUqY] ' 



and 



A(x, z, n) = + (z — uHq)^ 



(12) 



This solution to Laplace's equation, Eq. (8), does not, 
however, satisfy the boundary condition at the planar 
interface, Eq.(3), i.e. T^^-* ^ 1 at z = 0. As shown 
in 4], the isotherm at z = is deformed if the thermal 
conductance ratio a ^ 1. This, in turn, implies that the 
interface profile, which conforms to the isotherm at z = 0, 
is also deformed. However, in order to conduct a linear 
stability analysis of the planar interface, we must seek so- 
lutions to Eq. (7), with e = 0, and corresponding bound- 
ary conditions, Eqs.(3)-(6), that admit a planar interface 
profile. This is achieved by making use of the method of 
images. We refer the interested reader to ref. ^3 and 
references therein for more details on this method. We 
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FIG. 2: A sketch of a highly heat conducting inclusion in the 
melt (continuous line) and its image (dotted line), and the 
associated deformed solid-liquid interface profile. 



start by introducing an image inclusion that is centered 
at (0, —Hq). On using the principle of superposition, the 
equation 



T'j^\x,z) = 1 + Glz + 



(l-a)GL 
l+Co(a-l) 



(z - Ho)Cxi^x,z)\ 



[{z + Ha)Cx[x,z) + 



(13) 



not only satisfies Laplace's equation, but also the bound- 



ary condition, T 



1 at z = 0. This situation is de- 



picted schematically in figure 2 for the case of an inclusion 
and its image that are both characterized by a thermal 
conductance ratio a > 1. According to P, QjSH) the 
inclusion at (0, Hq) induces a depression in the interface 
shape while its image at (0, —Hq) induces a bump. The 
depression and the bump, being symmetrical about the 
a;— axis, cancel out leaving the interface planar. Note, 
however, that Eq. (13) fails to satisfy the boundary con- 
dition at the inclusion's surface dl, i.e. Eq.(5). Thus, It 
is necessary to place another fictitious particle at (0, 2Ho) 
and add its contribution to the thermal field. The result- 
ing equation will then satisfy Eq. (5) but fails to satisfy 
Eq. (3). So an image inclusion is placed at (0, —2Ho) and 
so on and so forth. This process is continued indefinitely 
to yield, 



(1 - a)GL 
l + Co(a-l) 



J2 iz + nHo)C, 



X{x.,z,n) • 



(14) 



n— — oo 



The expression, Eq.(14), now satisfies Laplace's equation 
and both boundary conditions, Eq.(3), at the planar in- 



terface, and Eq.(5) on dl. The thermal fields in the in- 
clusion and the solid phase are described by 



Gl{z~Ho) 
H-Co(a-l)' 



(15) 



T'j^\x,z) = 1 + Gsz + 



(1 - a)GL 
k[l + C^{a-\)] 



{z + nHa)Cx{x,z,n)- 



(16) 



n— — cc 
n#0 



(5) is) 

Note that Tg satisfies Laplace's equation, and = 1 
at the planar interface z = 0. The last term in Eq. (16) 
is necessary in order for Tg ' to satisfy the heat balance 
equation at the planar interface, i.e. the planar interface 

(S) 

moves with constant velocity; does not need to 

satisfy the boundary condition, Eq. (5), on the particle's 
surface dl . 



We now examine the linear stability of the planar state 
defined by Eqs. (14)-(16) and corresponding boundary 
conditions at the planar interface. We first superim- 
pose two-dimensional, time-dependent infinitesimal dis- 
turbances 6'(-^\ 0^^^ and rj upon the basic state solutions 



T; 



T)^' and z — 0, respectively in Eq. (7) (without the 



(L) rp{S) 

convection term) and boundary conditions, Eqs. (3)-(6). 
The resulting equations are then linearized with respect 
to the disturbances to yield the following, 



A6l(«) = 0, 



9'^^^ = -GlV + (^-q^ - F{x)r], at z = 0, (18) 



(17) 



6i(^) = -Gsv + a 



d^ri Fix) 
dx^ k 



rj, at z = 0, (19) 



dz 



0, as z ^ oo. 



d9^ 

dz 



0, as z ^ —oo, 
(20) 



-f^fc— J—, at z = 0. 21 

at dz dz 

The small slope approximation in the expression for the 
curvature, i.e., k « —d^rj/dx'^ has been used in Eqs. (18)- 
(19), and 



(1 - a)G, 
l+Co(a-l) ^ 



°^ dz +^^("-" 



(22) 
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The stability problem, Eqs. (17)-(21), is solved using the 
Fourier transform. We let T{f{x, z, t)} = f{uj, z, t). The 
transformed problem is 



52^(9) 



(23) 



with the accompanying conditions 
0<-^^ = -{GL + auj^)fi~T{F{x)Tj}, at z = 0, (24) 



k 

and with vanishing temperature gradients far away from 
the interface for ^'^^^ and 6'^"''\ In order to make progress, 
we consider the approximation T{F{x)ri} w F(0)ri since 
the major contribution to F{x) comes from near the ori- 
gin (see 'sl for a similar approximation for the case of a 
spherical particle), where F{0) can be be approximated 
by the leading asymptotic term for c « 1 (small deviation 
from a circular cylinder) as follows, 



FiO) 



{a — 1)Glc 1 
2[1 + Co{a - 1)]H^ ^ n 



2+0(c-l). (26) 



The approximation is, however, exact for a circular cylin- 
der (c = 1). On using the fact that X^^^i 1/"-^ = tt^/G 
and Co = l/(c + 1), we obtain 



F{Q) 



7r2(a-l)GLc(l + c) 



The solution of the transformed problem yields, 
e^^^ = -[Gl+ ou^^ + i^(0)]e-'^^ 



(27) 



(28) 



(29) 



On imposing the heat balance equation at the interface, 
we obtain the following evolution equation for the inter- 
face perturbation. 



di] 

at 



(30) 



where 17(w) = -{Q^{:\^k)aup- ^2F{<S)\uj, Q = kGs + GL 
is the conductivity weighted thermal gradient, and rep- 
resents the growth rate of the interfacial perturbation. If 
VL{lij) < then any initial disturbance dies out after a 
long time, and the planar interface is stable. If ^{llj) > 
then the opposite scenario holds, and the interface is un- 
stable. The marginal stability criterion is obtained by 
setting r2(w) = 0. A plot of the growth rate as a func- 
tion of the wavenumber lo is shown in figure 3 as func- 
tion of the aspect ratio, c, and the inclusion-interface 
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FIG. 3: A plot of the growth factor f2 versus the wavenumber 
uj for three values of the aspect ratio c (plot on left), and for 
three values of the inclusion-interface distance Ho (plot on 
right); ujc is the critical wavenumber, Eq. (31), and uin is the 
neutral mode, i.e. Q{ujn) ~ 0. The critical distance, He, is 
represented by a continuous line on the plot on the right. 



distance, Hq- Note that the range of unstable wavenum- 
bers, < UJ < LUn, increases with the aspect ratio c, where 
ujn is the neutral mode, i.e. 0(w„) = 0. Therefore, in- 
clusions with larger c values have a more destabilizing 
effect, or equivalently, the instability is onset at a larger 
inclusion-interface distance. The maximum growth rate 
occurs at the critical wavenumber ujc that is obtained by 
maximizing Q{uj) with respect to uj. We find. 



to. 



2F(0) + g 
3(H-fc)cr' 



(31) 



with the understanding that the parameter values are 
such that the expression in Eq. (31) is positive, i.e. 
F{0) < and (2F(0) + G) < 0. On substituting Eq. 
(31) into the expression for the growth rate, ri(w), we 
obtain the stability condition, Q < — 2F(0), and on using 
Equation (27) , we find the following stability criterion for 
an elliptic cylinder, valid for c w 1, 



ln^c{l+c)GL{l-a) 



(32) 



y 6ic + a)g 

which reduces to the following exact criterion for a cir- 
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cular cylinder, 
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(* 1.0e-04) 



FIG. 4: Plot of the critical inclusion-interface distance, He 
(Eq. (32)), as function of the thermal gradient Gl (note the 
scale factor l.Oe — 04) for three distinct values of the aspect 
ratio c, c — 1.2, c = 1.0 and c = 0.8, for a system consisting 
of succinonitrile (SCN) containing an inclusion with (a = 5) . 
The regions of instability in the Hq — Gl plane are below the 
curves. The numerical values of the parameters used in the 
calculations are taken from Ref. [Tl|. 



He 



In-^GLil-a) 
3{l + a)g ■ 



(33) 



The critical distance He is the largest value of Hq for 
the instability to appear. Note that the stability con- 
ditions, Eqs. (32)-(33), require that (i) the product 
G'l(1 - a) > 0, and (ii) He > 1. Therefore, the in- 
stability is observable only for certain combinations of 
the physical and processing parameters. Figure 4 depicts 
the plot of He versus the dimensionless thermal gradient 
in the liquid for three distinct values of c. The numerical 
values of the parameters used in calculating He pertain 
to an experimental system consisting of SCN containing 
an inclusion whose thermal conductance ratio, a, is ar- 
bitrarily set equal to 0.5. We note the following: (i) He 
increases with the aspect ratio c, (ii) He increases with 
Gl, and (iii) for every value of the aspect ratio c, there 
exists a lower bound for Gl below which the instability 
is not observable, i.e. He < 1. The dependence of He on 
the inclusions size, i.e. a, and on other factors, such as a, 
resemble that of the case of a spherical particle, with the 
exception that for spherical inclusions. He ^ (Gl/GY^^- 
We refer the interested reader to 6] for more details. 
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